ASTROPHYSICAL JOURNAL LETTERS, SUBMITTED 

Preprint typeset using LfT^C style emulateapj v. 14/09/00 



O 

o 



> 
in 

in 
in 
o 



^ ■ 

Or 

6 ■ 

a: 



13 



THE SUPER STAR CLUSTER NGC 1569-A RESOLVED ON SUB-PARSEC SCALES 
WITH HUBBLE SPACE TELESCOPE SPECTROSCOPY 1 

Dan Maoz 2 3 , Luis C. Ho 4 , and Amiel Sternberg 3 

Received 2001 May 1; Accepted May 22 

ABSTRACT 

We present 3000-10000 A HST/STIS long-slit spectroscopy of the bright super star cluster A (SSC-A) in the 
dwarf starburst galaxy NGC 1569. The 0"05 HST angular resolution allows, for the first time, to probe for spatial 
variations in the stellar population of a ~ 10 M SSC. Integrated ground-based spectra of SSC-A have previously 
revealed young Wolf-Rayet (WR) signatures that coexist with features from supposedly older, red supergiant 
(RSG), populations. We find that the WR emission complexes come solely from the subcluster A2, identified 
in previous HST imaging, and are absent from the main cluster Al, thus resolving the question of whether the 
WR and RSG features arise in a single or distinct clusters. The equivalent widths of the WR features in A2 — 
including the C IV A5808 complex which we detect in this object for the first time — are larger than previously 
observed in other WR galaxies. Models with sub-solar metallicity, as inferred from the nebular emission lines of 
this galaxy, predict much lower equivalent widths. On the "clean" side of Al, opposite to A2, we find no evidence 
for radial gradients in the observed stellar population at 0"05< R <0."40 (~ 0.5 to 5 pc), neither in broad-band, 
low-resolution, spectra nor in medium-resolution spectra of the infrared Ca II triplet. 

Subject headings: galaxies: individual (NGC 1569) — galaxies: star clusters — galaxies: starburst — galaxies: 
stellar content — stars: Wolf-Rayet 



1. INTRODUCTION 

A substantial fraction of the star formation in nearby star- 
bursts, be it in merging galaxies, dwarf galaxies, or early-type 
galaxies with circumnuclear rings, takes place in so-called su- 
per star clusters (SSCs; e.g., Whitmore et al. 1999; Hunter et 
al. 2000; Maoz et al. 2001, and references therein). These 
are clusters of stars with total luminosities as high as My = -15 
mag (Ly = 1.3 x 10 8 Lyo) and radii of order a few pc. SSCs 
present a mode of massive star formation distinct from the OB 
associations found in most star-forming regions of quiescent 
galaxies. They are of interest for understanding star formation 
at high redshift, which is dominated by the starburst mode, and 
the formation of globular clusters, into which SSCs may evolve 
if they contain sufficient numbers of low-mass stars, and can 
survive disruption and evaporation. 

Two of the nearest (2.5 ±0.5 Mpc; O'Connell, Gallagher, 
& Hunter 1994) and best-studied SSCs are clusters A and B 
in the amorphous dwarf galaxy NGC 1569 (Arp & Sandage 
1985; Melnick, Moles, & Terlevich 1985). NGC 1569-A is one 
of only several SSCs that have a kinematic mass estimate, w 
1O 6 M (Ho & Filippenko 1996; Sternberg 1998). Prada, Greve, 
& McKeith (1994) detected the Ca II AA8498, 8542, 8662 ab- 
sorption triplet, which forms in the atmospheres of red giants 
and supergiants, in cluster A. Based on the equivalent width 
(EW) of these lines, they estimated an age > 13 Myr. Gonzalez 
Delgado et al. (1997), however, detected in the SSC-A spec- 
trum the bump around 4686 A generally ascribed to He II, C III, 
C IV, N III, and N V emission from Wolf-Rayet (WR) stars. WR 
stars are the evolved stages of stars more massive than ~ 30M Q , 
and are thought to exist only when a burst is 3 - 5 Myr old. 
To explain this, Gonzalez Delgado et al. (1997) proposed that 
there have been several succesive bursts of star formation in 



SSC-A. De Marchi et al. (1997) analyzed post-refurbishment 
Hubble Space Telescope (HST) images and showed that SSC-A 
has two peaks in its light distribution, separated by 0" 1 8, which 
they designated Al and A2, and postulated that these are two 
distinct clusters, with the WR features coming from one cluster 
and the red supergiant signatures from the other. More recently, 
Hunter et al. (2000) obtained deeper HST images in several 
bands, and argued against the proposed solution of De Marchi 
et al. (1997), based on the finding that Al and A2 have sim- 
ilar colors. They noted that WR stars and red supergiants are 
expected to co-exist, to some degree, when a cluster is about 
5 Myr old. Nonetheless, Buckalew et al. (2000) have found ev- 
idence, based on HST He II narrow -band imaging, that the WR 
emission is concentrated in A2. 

More generally, how SSCs form and survive as bound ob- 
jects is an unsolved problem, given the rapid ejection of gas 
from the proto-clusters by jets, ionizing radiation, stellar winds, 
and subsequent supernovae. Some clues may be found in the 
spatial structure of the nearest SSCs. NGC 3603 in the Milky 
Way and R136 in the Large Magellanic Cloud (LMC) show 
some evidence for segregation of the stellar population, with 
the most massive stars concentrated at smaller radii (Moffat, 
Drissen, & Shara 1994; Brandl et al. 1996; Hunter et al. 1996). 
It is thought that the clusters are too young for dynamical mass 
segregation to have occurred (Bonnell & Davies 1998). The 
segregation therefore may reflect a gradient in age or in ini- 
tial mass function (IMF) properties, but could also be a sta- 
tistical effect, due to the rarity of high-mass stars in the outer, 
low stellar-density regions of the clusters. Current estimates for 
the masses of these clusters are ~ 10 4 - 10 5 M Q (Drissen et al. 
1995). The spatial population profile of massive, 10 6 M Q SSCs 
such as NGC 1569-A has never been observed. This Letter 



presents resolved HST long-slit spectroscopy of NGC 1569-A 
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to address several of these issues. 

2. OBSERVATIONS AND REDUCTION 

We observed NGC 1569 on 28 November 2000 UT, using 
the Space Telescope Imaging Spectrograph (STIS) on HST . 
The 52" x 0."1 slit was centered on SSC-A, and the spacecraft 
was oriented with the intention that clusters Al and A2, and 
cluster B, 7" to the south-east, would all fall within the narrow 
slit. Unfortunately, a 0.°45 error in the header orientation infor- 
mation of a WFPC2 image of the galaxy from 1996, combined 
with a 0.°33 error in determination of the STIS slit orientation 
in the focal plane at the time the observation was planned, re- 
sulted in a total orientation error of 0.°78. The 0"l-wide slit 
was misplaced 0"1 southwest from the center of SSC-B, effec- 
tively missing it. A repeat observation of SSC-B is planned, 
and in this Letter we will analyze only the spectra of SSC-A. 




Fig. 1- Section (l"15x If 15) of STIS acquisition image of 
NGC 1569-A. The slit position is shown, and the various spa- 
tial bins for which spectra were extracted are labeled. 



Exposure times totaled 32.4 min in the G750L grating 
(5250-10250 A, 10 A resolution full width at half maximum 
[FWHM]), 43.2 min in the G430L grating (2910-5714 A, 5.5 
A resolution), and 93.6 min in the G750M grating (8276-8842 
A, 1.1 A resolution). The STIS plate scale is 0"05 pixel -1 . 
Flat-field exposures were taken with the G750L and G750M 
gratings for use in fringe correction. The data were initially 
processed by the STScI pipeline, producing geometrically cor- 
rected, flux- and wavelength-calibrated two-dimensional spec- 
tral images at each of three dither positions. The fringe flats 
were normalized and then scaled to give the best suppression 
of the fringing in the science images. The images at the indi- 
vidual dither positions were registered and combined. Spectra 
of each of the resolved spatial bins of SSC-A were then ob- 
tained by summing individual rows. The five spatial bins cho- 
sen were a central bin, consisting of the two rows with peak 
flux (R <0."05), two intermediate-radius bins, each consisting 
of two rows above and below the central bin (0."05< R <0"15), 
and two outer bins (0."15< R <0."40). The outer bin on the 
south-east side of SSC-A is dominated by flux from SSC-A2, 
while the center and north-west bins are dominated by flux from 
SSC-A1. Figure 1 shows the location of the rectangular sec- 
tions of SSC-A we have thus probed. 



Despite the geometric correction applied to the data in the 
pipeline, some residual curvature remains in the spectral im- 
ages. Combined with the unavoidable variations in point-spread 
function (PSF) along the large wavelength interval of the low- 
resolution spectra, this produces conspicuous artificial spectral 
variations between different rows. To correct for these effects, 
we obtained the G430L and G750L archival data of two white 
dwarfs, GRW+70D5824 and Feige 110, respectively, which 
were observed with the same STIS setup. We simulated an 
A1+A2 cluster pair having no spatial color gradient by con- 
volving the spectral images of the stars with Gaussians, scal- 
ing, shifting, and adding the images, according to the widths 
and relative fluxes measured in WFPC2 images by Hunter et 
al. (2000). We verified that the cross-dispersion profile at 
5500 A of the simulated data matched closely that of the real 
data, which have FWHM of 0"20±0"02. We then extracted 
spectra in the same spatial bins as above, normalized by the to- 
tal spectrum of the simulated cluster. The resulting correction 
function was heavily smoothed to eliminate residual spectral 
features from the white dwarfs, and each spatial bin in the real 
data was divided by its corresponding correction function. The 
correction functions are insensitive, at a level of a few percent 
or less, to the details of the simulation, as long as it includes 
2 sources with approximately the correct scaling and the total 
profile has 0."15<FWHM<0"27. 

Background emission from the sky and from the galaxy in the 
neighborhood of the cluster is negligible. To obtain a "clean" 
spectrum of A2, the spectrum of the Al bin diametrically op- 
posed to it was subtracted from the A2 bin. Residual cosmic 
rays and bad pixels in the final spectra were identified by eye 
and corrected by interpolation between adjacent pixels. There 
is agreement in slope and flux level between the G430L and 
G750L spectra in the overlap region at 5250 A to 5714 A to 
the few-percent level, except for the central bin spectra, where 
the G430L spectrum must be scaled up by 15% to match the 
G750L data. After accounting for the 40% transmission of the 
total light by the 0."1 slit, the total SSC-A flux matches well 
the measurements by Gonzalez Delgado et al. (1997) and De 
Marchietal. (1997). 

3. ANALYSIS AND DISCUSSION 
3.1. Subcluster A2 and the Wolf-Rayet Emission 

Figure 2 shows the final low-resolution spectra for each of the 
spatial bins. To allow easy comparison of the spectral shapes, 
the spectra have been multiplicatively scaled up to match the 
central-bin spectrum in the blue region, and then shifted up or 
down by an additive constant. The spectrum of the A2 bin 
is shown before subtraction of the northwest, diametrically- 
opposed, bin of Al (NW 0."15< R <0."40). The EWs of the 
main emission and absorption features are listed in Table 1 . 

Gonzalez Delgado et al. (1997) first reported the WR emis- 
sion bump at 4686 A in the integrated spectrum of SSC-A. It 
is clear from Figure 2 that this feature, as well as the strong 
C IV A5808 A WR complex which we identify for the first time 
in this object, come almost exclusively from SSC-A2. Indeed, 
the weakness or absence of the WR features in the Al bins 
demonstrates the angular resolving power of this observation. 
These results resolve the issue of the co-existence, in SSC-A, 
of WR with red supergiant signatures, confirming the interpre- 
tation proposed by De Marchi et al. (1997) and the narrow- 
band imaging results of Buckalew et al. (2000). It is unknown 
whether, at the 0" 1 8 (2.2 pc) separation, A2 is a separate cluster 
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seen in projection, or rather a physical neighbor or substructure 
to Al. In either case, the data show that A2 constitutes a dis- 
tinct stellar population. Note that, apart from the WR features 
in A2, the stronger Balmer-line emission in A2, and a possible 
difference in the spectral shape at A < 3600 A, the spectra of 
Al and A2 are quite similar. This confirms and explains the 
observation by Hunter et al. (2000) of similar optical colors in 
Al and A2. 
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Fig. 2- Broad-band spectra for each of the spatial bins of 
NGC 1569-A. The flux scale corresponds to the second spec- 
trum from the top (the central, R <0"05, bin of Al). Note 
the strong WR emission that comes almost exclusively from 
the A2 bin, and the similarity of the continuum shape among 
all bins. Scaling and shifting that has been applied to the 
spectra, from top to bottom, is 2.5/ A + 2.2 X 10~ 16 ; fx +0; 
2.25A-3.5 x 1CT 16 ; and 5.6/ A -6.5 x lCr 16 . 

The flux in the 4686 A feature is consistent with the measure- 
ment by Gonzalez Delgado et al. (1997), who concluded it re- 
quires the luminosity of 20-40 WNL stars. However, now that 
this emission is localized to the relatively faint sub-clump A2, 
the EW of the WR features becomes huge. After subtraction of 
the background from Al, we measure EW(4686)=33±9 A and 
EW(5808)=32±9 A. For comparison, Schaerer, Contini, & 
Kunth (1999) have measured in five WR galaxies typical val- 
ues of EW(4686)=4-10 A and EW(5808)=2-6 A. They find the 
highest values, with EW(5808)~ 12 A, in the galaxy Tol 89. 
Schaerer & Vacca (1998) have constructed evolutionary syn- 
thesis models for instantaneous starbursts and predicted the 
strengths of the WR features as a function of starburst age and 
metallicity. As seen, for example, in Figure 1 1 of Schaerer 
& Vacca (1998), the WR emission peaks briefly between ages 
of 3-5 Myr, and its strength is proportional to the metallicity 
Z. For Z = 0.2Z Q , they predict that the EW of either of the 
WR bumps never rises above 8 A, and for Z = 0.4Z Q , it has 
a maximum of 1 1 A. Only at Z = 2Z Q do the EWs ever reach 
EW(4686)=30 A and EW(5808)=25 A. 

However, several measurements of the gas chemical abun- 
dances in NGC 1569 have yielded a significantly sub-solar 
metallicity, of Z w 0.25Z Q (Calzetti, Kinney, & Storchi- 
Bergmann 1994; Kobulnicky & Skillman 1997; Gonzalez Del- 
gado et al. 1997). The large EWs of the WR features we 



find in A2 would therefore indicate that, contrary to the line- 
emitting gas from which they were recently formed, the stars 
in A2 have solar or supersolar abundances. Alternatively, there 
may be a problem with the models, in that the WR lifetimes 
depend on uncertain mass-loss rates and the effects of stel- 
lar rotation (Maeder & Meynet 2000). Note that, in a non- 
instantaneous starburst, the EW of the WR features will nec- 
essarily be smaller, thus exacerbating the problem. If A2 is not 
a separate cluster, or a subcluster of A, but rather a particular 
physical region of cluster A, then it is not clear that subtraction 
of the Al background is justified. In that case, the EWs (and 
the errors) of the WR features are halved, but are still double 
the values expected for the metallicity of the gas. 

3.2. Cluster Al Population and Radial Color Gradient 

The data for the "clean," northwest, side of SSC-A1 in Fig. 
1 reveal the spectrum of Al, uncontaminated by A2. We de- 
fer a detailed analysis and comparison with models to a future 
paper. However, except for the WR features which we have 
shown come from A2, the main characteristics of the spectrum 
are unchanged in comparison to the integrated spectra analyzed 
by previous ground-based studies (e.g., Gonzalez Delgado et 
al. 1997). Specifically, the blue slope and the weakness of the 
Balmer jump still point to an age of ~ 5 Myr, and the presence 
of the Ca II triplet absorption rules out a younger age. 

It is remarkable that Al and A2 appear to have similar ages 
of ~ 5 Myr, but that one has WR emission and the other does 
not. A possible explanation is that Al has an IMF with an up- 
per cutoff <, 30M Q , and therefore does not form WR stars. A 
higher mass cutoff in A2 could be related to its anomalously 
high metallicity, suggested by the large EWs of its WR fea- 
tures. Alternatively, Al may also have high, approximately 
solar, metallicity, but is slightly older than A2, say ~ 7 Myr, 
by which time it would have exited its WR phase (Schaerer 
& Vacca 1998). Finally, Al may have the normal, sub-solar 
metallicity of the galaxy, in which case it would just be exiting 
its WR phase already at 5 Myr. Thus Al and A2 must have 
widely different IMFs, or widely different abundances, or sim- 
ilar, anomalously high, abundances but slightly different ages. 

The northwest Al data also allow us to examine, for the first 
time, the radial gradient in the stellar population of an individ- 
ual high-mass SSC. As seen in Figure 2, there is remarkable 
similarity in the spectra of Al in the three radial bins shown. 
This indicates that, at least on its northwest side, Al is ex- 
tremely homogeneous. Considering the small age, the forma- 
tion of stars in Al must have been highly synchronized through- 
out its volume. In this respect, SSC-A1 appears to differ from 
other, lower mass, SSCs that possibly do show population seg- 
regation (see §1). 

3.3. Calcium Triplet Spectra 

Figure 3 shows a region of the medium-resolution G750M 
spectra centered on the near-infrared Ca II triplet. The spec- 
tra of the outer bins are too noisy to reveal unambiguously the 
Ca II triplet, and are not shown. Figure 3 and Table 1 show 
that there is no significant difference in Ca II triplet absorption 
between the central and intermediate radial bins in SSC-A1. 
This result confirms the lack of a population gradient seen in 
the low-resolution, broad-band spectra. Although the strength 
of the Ca II triplet absorption, which comes first from red su- 
pergiants and at a later stage from red giants, has been used to 
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argue for a large age for SSC-A (Prada et al. 1994), the data are 
also consistent with an age of ~ 5 Myr. Garcia- Vargas, Mollan, 
& Bressan (1998) have used evolutionary synthesis models to 
predict the EW for the sum of the two strongest triplet lines, at 
8542 A and 8662 A. The summed EW value seen in NGC 1569, 
of ~3-5 A, is generally lower than the value predicted for star- 
bursts older than 6 Myr, regardless of metallicity (cf. Boker 
et al. 2001 for caveats on the accuracy of the low-metallicity 
models). The observed EW, on the other hand, briefly matches 
the model predictions for cluster ages 5-6 Myr, when the red 
supergiants first form. 
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Fig. 3- Medium-resolution spectra centered on the near- 
infrared Ca II triplet absorption. The lower spectrum, which 
is the sum of the spectra of the two bins straddling the central 
bin, has been shifted down by 0.8 units. Note the similarity of 
Ca II triplet absorption strength at the different radii. 



3.4. Summary 

We have obtained spectra of NGC 1569- A that are resolved 
at the sub-pc scale, the first time a 10 6 M Q SSC is observed this 
way. The data clearly show that the WR emission from this 
SSC is confined to subcluster A2. The EWs of the WR features 
are much larger than previously observed in other WR galaxies, 
and are several times larger than the maximum values predicted 
from models with sub-solar metallicities, as measured in the gas 
of this galaxy. This problem is particularly severe if A2 is a sep- 
arate structure (as opposed to a region inside Al), justifying the 
subtraction of the underlying light of Al from A2. The con- 
tinuum slope, the Balmer jump, the Ca II absorption and WR 
emission EWs, are all consistent with an age of ~ 5 Myr for 
the bursts that formed both Al and A2. While there is a sharp 
difference in WR population between Al and A2, there is no 
significant spatial gradient in the stellar population on the NW 
side of Al. 

We conclude that, either A2 has an anomalous, super-solar 
metal abundance, or the models have seriously underpredicted 
the WR emission strength for sub-solar metallicities. The 
weakness of WR features in Al, despite the similar age to A2, 
could be caused by a much lower IMF upper-mass cutoff in Al, 
a much lower metallicity in A 1, or a slightly larger age for Al. 
Finally, the homogeneity of Al is distinct from the possible ra- 
dial mass segregation reported in some lower mass SSCs in the 
Galaxy and the LMC. Further analysis is required to determine 
the degree to which the apparent uniformity may be set by the 
most massive stars, which dominate the light. If real, it needs 
to be reproduced by models of SSC formation and evolution. 

This work was supported by grant GO-8293 from the Space 
Telescope Science Institute, which is operated by AURA, Inc., 
under NASA contract NAS 5-26555. 



Table 1 

ngc 1569- a: emission and absorption equivalent widths 



A2 


Al 






R<0"05 0."05<fl <0."15 


0."15<fl <0"40 



H5A4102 


-0.5 ±0.5 


-2.8 ±0.4 


-3.0 ±0.7 


-1±1 


H 7 A4340 


1.2±0.8 


-2.0 ±0.3 


-2.5 ±0.5 


-2±1 


WR A4686 


33±9 


1.9 ±0.5 


1.1 ±0.5 


< 1.5 


H/3A4861 


11±1 


-0.3 ±0.1 


-0.9 ±0.3 


0±1 


WR A5808 


32±9 


<3 


<3 


<8 


HaA6563 


56±6 


8.8±0.5 


8.5 ±1.0 


10±2 


Ca II A8498 




-1.1 ±0.2 


-0.9 ±0.2 




Ca II A8542 




-2.4 ±0.6 


-1.5±0.5 




Ca II A8662 




-1.7 ±0.5 


-1.7 ±0.6 





Note. — Equivalent widths in A. Positive values are for emission and negative 
for absorption. WR denotes the broad Wolf-Rayet emission complexes roughly 
centered on the wavelengths listed. Errors correspond to ~ 2er and are dominated 
by continuum placement uncertainty. 



5 



REFERENCES 



Arp, H., & Sandage, A. 1985, AJ, 90, 1163 

Boker, T., van der Marel, R. P., Mazzuca, L., Rix, H.-W., Rudnick, G., Ho, 

L. C, & Shields, J. C. 2001, AJ, 121, 1473 
Bonnell, I. A., & Davies, M.B. 1998, MNRAS, 295, 691 
Brandl, B., et al. 1996, ApJ, 466, 254 

Buckalew, B. A., Dufour, R. J., Shopbell, P. L., & Walter, D. K. 2000, AJ, 120, 
2402 

Calzetti, D., Kinney, A. L. & Storchi-Bergmann, T. 1994, ApJ, 429, 582 

De Marchi, G., Clampin, M., Greggio, L., Leitherer, C., Nota, A., & Tosi, M. 

1997, ApJ, 479, L27 
Drissen, L., Moffat, A. F. J., Walborn, N. R., & Shara, M. M. 1995, AJ, 110, 

2235 

Garcia- Vargas, M. L., Molla, M., & Bressan, A. 1998, A&AS, 130, 513 
Gonzalez Delgado, R. M., Leitherer, C, Heckman, T., & Cervifio, M. 1997, 

ApJ, 483, 705 
Ho, L. C, & Filippenko, A. V. 1996, ApJ, 466, L83 
Hunter, D. A., O'Connell, R. W., & Gallagher, J. S. 1994, AJ, 108, 84 



Hunter, D. A., O'Connell, R. W., Gallagher, J. S. & Smecker-Hane, T. A. 2000, 
AJ, 120, 2383 

Hunter, D. A., O'Neil, E. J., Lynds, R., Shaya, E. J., Groth, E. J., & Holtzman, 

J. A. 1996, ApJ, 459, L27 
Kobulnicky, H. A. & Skillman, E. D. 1997, ApJ, 489, 636 
Maeder, A. & Meynet G. 2000, ARAA 38, 143 

Maoz, D., Barth, A. J., Ho, L. C, Sternberg, A., & Filippenko, A. V. 2001, AJ, 

in press ((astrop-ph/0103213) 
Melnick, J., Moles, M., & Terlevich, R. 1985, A&A, 149, L24 
Moffat, A. F. J., Drissen, L., & Shara, M. M. 1994, ApJ, 436, 183 
O'Connell, R. W., Gallagher, J. S., & Hunter, D. A. 1994, ApJ, 433, 65 
Prada, F, Greve, A., & McKeith, C. D. 1994, A&A, 288, 396 
Schaerer, D., Contini, T., & Kunth, D. 1999, A&A, 341, 399 
Schaerer, D. & Vacca, W. D. 1998, ApJ, 497, 618 
Sternberg, A. 1998, ApJ, 506, 721 

Whitmore, B. C, Zhang, Q., Leitherer, C, Fall, S. M., Schweizer, F. & Miller, 
B. W. 1999, AJ, 118, 1551 



